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In layered single band systems, the interlayer conductivity depends on the orientation of 
the in-plane magnetic field and takes maximum values when the magnetic field is perpendic- 
ular to flat regions of the Fermi surface. Extending this known results to multi-band systems, 
we propose an experiment to extract information about their Fermi surface topology. We dis- 
cuss application of the formula to a FeAs-based superconductor, LaFeAsOi_a;Fj;. We show 
that the magnetically ordered state in the parent compound is clearly distinguished from 
the paramagnetic state by the oscillation period in the interlayer conductivity. We demon- 
strate that evolution of the Fermi surface topology by changing the doping concentration is 
reflected to the interlayer conductivity oscillation patterns. 
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1. Introduction 

Determination of the Fermi surface topology plays an important role to make clear the basic elec- 
tronic structure of a superconductor. For layered compounds, angle-resolved photoemission spectroscopy 
(ARPES) can be used to study the Fermi surface. Although ARPES provides the map of the Fermi surface, 
ARPES measurements are sensitive to surface conditions. Contrary, magnetic quantum oscillations provide 
information about the Fermi surface area perpendicular to the applied magnetic field. The Fermi surface 
areas determined by magnetic quantum oscillations are bulk electronic properties but this technique does 
not provide information about the Fermi surface topology directly. In this paper, we propose interlayer 
magnetoresistance experiments to probe the bulk Fermi surface topology in layered multi-band systems 
and its application to the newly discovered FeAs-based layered superconductors.^ To be specific, we focus 
on LaFeAsOi-^^F^: and discuss what we expect for the interlayer magnetoresistance in the paramagnetic 
state and the magnetically ordered state. 

The Fermi surfaces of the FeAs-based layered superconductors have been studied theoretically and 
experimentally. First principles calculations^'^ suggest that the Fermi surface of LaFeAsO consists of hole 
Fermi surfaces around F point and electron Fermi surfaces around M point. The Fermi surfaces measured 
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by ARPES*^"^ are mostly consistent with this picture, though the measured samples are limited. For 
example, a single hole pocket and a single electron pocket were observed in NdFeAsOo.gFo.i-j^''^ while 
the other Fermi pockets predicted by the band calculations were not observed. The hole Fermi surfaces 
and the electron Fermi surfaces have instability of a magnetic order associated with a nesting of (vr, vr) 
wavevector. The magnetically ordered state is a collinear antiferromagnetic state as suggested from the 
neutron scattering experiments."'^''' ^-"^ Quantum oscillation measurements were carried out for magnetic 
states of a few Fe-based superconductors. ^^"^^ Possible folding of Brillouin zone are discussed so that the 
observed oscillation periods are reproduced. 

Another experiment to probe the Fermi surface of layered compounds is to measure angular mag- 
netoresistance oscillations which was firstly discovered in the organic conductors /3-(BEDT-TTF)2lBr2"'^^ 
and 0-(BEDT-TTF)2l3.^^ The magnetoresistance oscillations in these systems were well understood by 
a semiclassical analysis assuming an energy dispersion described by a two-dimensional cylinder with a 
weak cosine warping.^''' In this magnetoresistance phenomenon, if the applied magnetic field lies in the 
plane, ^^"^^ the oscillation patterns are associated with the Fermi surface topology. The interlayer conduc- 
tivity is dominated by the regions of the Fermi surface where the Fermi velocity is parallel to the magnetic 
field. The theoretical calculation for Tl2Ba2Cu06^^ is consistent with the experiment. 

In this paper, we extend the formula for the single-band systems^^"^"*^ to multi-band systems. The 
derivation is given in Sec. 2. We apply the formula to the paramagnetic state and the magnetic state of 
LaFeAsO (Sec. 3) and the paramagnetic state of LaFeAsOi-ajF^: (Sec. 4). The electronic band structure of 
the parent compound is determined by first principles calculations. For doped systems, we use the five 
band model described in ref.23. We consider LaFeAsOi_a;Fa; because two-dimensionality is stronger than 
other FeAs-based superconductors, such as BaFe2As2. 

2. Interlayer Conductivity Formula 

In this section, we derive the angular magnetoresistance formula for multi-band systems. We assume 
that the system has a structure of stacking conduction layers in the z axis. The magnetic field B = 
{Bx,By,0) = B{cos(f),sm(l),0) is assumed to be in the plane. The Fermi surface topology is associated 
with the azimuthal angle (f) dependence of the interlayer magnetoresistance as we shall see below. The 
vector potential A is taken as A = {ByZ, —B^z, 0). 

The interlayer tunneling between j-th. plane and j + 1-th plane is described by 

Ht = -tc Yl I d'r^^l. (r) i^^^l (r) exp {ix (r)) + h.c, (1) 

where v represents the band index. Here we have included the effect of the in-plane magnetic field via a 
Pierls procedure: 

X (r) = Jrflc {Byx - B^y) , (2) 
with ttc the interlayer spacing along the z axis. The interlayer conductivity is calculated by the Kubo 
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formula, 

ImK^. (to) 

where K^^ (w) is the Fourier transform of the following retarded current-current correlation function, 

- ^0 = (* - [-^^ ' "^^ (*')]> ' (4) 

with S the area of the system. We calculate K^^{t) in the Matsubara formalism. The correlation function 
in the Matsubara formalism is given by 

K^z (r) = ^ {TrJ. (r) J. (0)) (5) 

Here the current operator is 

J. = '^Y.Jd'r (r) 41 (r) exp {ix (r)) + h.c. (6) 

The function K^^ (a;) is obtained by K^^ (w) = [iun oj + i6) where K^^ {ioJn) = 

dr exp (iLOnT) (r) and 5 is a positive infinitesimal number. For the non-interacting case, we ob- 
tain 

2Trh 



k,k',i/,j,(T 



x\{k',j + l,a\j,\k,j,a)\\ (7) 

where the i^-th energy band is denoted by Substituting the matrix element of the current operaotr 

Jz for the electron hopping from the j-th layer to the j + 1-th layer, 

(k',j + l,o-| |k, j,c7) = — ^5k,k'-b, (8) 
with b = eac{—By,Bx)/{ch), we obtain 

x^(e-eI:1)--^ . (9) 



Here we have replaced the second Dirac delta function in eq.(7) by a Lorentz function with F = ^/(2t) to 
include the scattering effect upon hopping between layers. Although it is possible to take different F's for 
each band, we take the same F for all bands for simplicity. The formula (9) is understood in a following 
way. According to the semiclassical equations under a magnetic field, an electron moves in k-space upon 
hopping between one layer to the adjacent layers. The contribution to a^z is large if the motion in k-space is 
along the Fermi surface because the Lorentz factor in eq. (9) takes the maximum value. Since the direction 
of the motion is perpendicular to the in-plane magnetic field, dominant contributions arise if the magnetic 
field is perpendicular to the fiat regions of the Fermi surface. 
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At zero temperature, the integration over the BriUouin zone is transformed to the integration along 
the Fermi surfaces as fohows 



where = {l/h)dE^ /d\i. From the form of the Lorentz factor, it is clear that dominant contributions 
to Gzz arise from the Fermi surface parts with the Fermi velocity being parallel to the magnetic field, which 
is consistent with the intuitive picture mentioned above. The weight of the contribution is determined by 
the factor l/lv^'^'*! as in the density of states. Anisotropy in (7 zz detects deviation of the Fermi surface 
shape from the circular shape. If a Fermi surface is completely circular, such a Fermi surface contributes 
to azz only as the constant value independent of the orientation of the magnetic field. In order to calculate 
azz from cq.(lO) wc need values of v^*^^ and the Fermi surface data for each band. For Fermi surfaces which 
are not given by separated lines, one needs to calculate cq.(9). 

Note that wc have neglected the Zeeman energy contribution in eq.(lO). To be precise, wc need the 
energy contour at i^^^) = Ep it gfisB. However, the Zeeman energy is much less than the Fermi energy. 
Therefore, we may neglect the Fermi surface change associated with the Fermi energy shift by the Zeeman 
energy. 

3. Interlayer Conductivity of LaFeAsO 

For the application to LaFeAsO, we compute the electronic band structure by first principles calcu- 
lations. The Fermi surfaces were obtained by the first-principles density functional calculation package 
VASP.^^ Figure 1(a) shows the Fermi surface of the paramagnetic state and Fig. 1(b) shows the interlayer 
conductivity for F = 2 x 10"^ and B = lOT which corresponds to uJcT = eBT/{mc) = 0.29. (Hereafter 
we measure F in units of eV.) Each band contribution to the interlayer conductivity is also shown. The 
values are normalized by the total interlayer conductivity at = 0°. The broad peaks in the interlayer 
conductivity at (f) = 45° and (f) = 135° are associated with the contributions from the Fermi surfaces FS3 
and FS5 because those Fermi surfaces have flat regions whose Fermi velocities are in these directions. 
The peak at cf) = 90° is associated with the flat region in the Fermi surface FS4. The Fermi surfaces FSl 
and FS2 are nearly circular. Therefore, their interlayer conductivities are almost featureless. The largest 
contribution to azz comes from FS2 but other band features are clearly seen because anisotropy in the 
FS2 contribution is small. 

The Fermi surface and the interlayer conductivity in the magnetically ordered state is shown in Figs. 2(a) 
and 2(b). The magnetic ordered state is a collinear antiferromagnetic order.^ The interlayer conductivity 
is dominated by the contribution from the Fermi surface around the F point (FSl). The broad peak at 
(f) = 90° is associated with the fiat regions in FSl. The small Fermi surfaces FS2 and FS3 are almost 
featureless because they are nearly circular. 
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Fig. 1. (Color online) The paramagnetic state Fermi surface (a) and the normalized interlayer conductivity (b). 
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Fig. 2. (Color online) The Fermi surface of the antifcrromagnetic state (a) and the normalized interlayer conduc- 
tivity (b). 



A clear difference between the paramagnetic state and the antiferromagnetic state is the oscillation 
period difference in azz- In the paramagnetic state, the oscillation period is 90°. By contrast the oscil- 
lation period of the antiferromagnetic state is 180°. Therefore, the presence of the magnetic ordering is 
distinguished from the paramagnetic state by measuring the oscillation period of Ozz- 

In order to make clear the condition of observing the oscillations in the interlayer conductivity, we 
compute the magnetic field dependence of the peak height for several values of F. The onset of the 
anisotropy depends on F. Therefore, the measurements of the magnetic field dependence of the anisotropy 
provides information about the scattering. Figure 3(a) shows the magnetic field dependence of the ratio of 
the peak in Uzz at 45° and the value at 0° in the paramagnetic state. The ratio increases monotonically as 
we increase the magnetic field. Figure 3(b) shows the magnetic field dependence of the ratio of the peak 
at 90° and the dip at 75.7°. Since the magnetic field effect on those angles are not simple, the ratio shows 
a non-monotonic behavior. In principle one can deduce F by measuring these two ratios by comparing 
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the theoretical prediction based on first principles calculations with the experiments. Figure 4 shows the 
magnetic field dependence of the ratio of the peak in azz at 90° and the value at 0° in the magnetically 
ordered state. 




(a) (b) 

Fig. 3. (Color online) Magnetic field dependence of a zz{^b°) / a zz{^° ){&) and (Jzz{^^°) / cfzz[^^■^°)iy) in the param- 
agnetic state. 
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Fig. 4. (Color online) Magnetic field dependence of <J zz{^^°) / zz{^°) in the magnetic state. 



4. Doped Systems 

In this section we consider the doping dependence of the interlayer conductivity. Our consideration is 
based on the rigid band picture, and we take the five-band model Hamiltonian from Kuroki et al?^ Figure 
5 shows the electron density dependence of the Fermi surface and Fig. 6 shows the interlayer conductivity 
for each electron concentration with B = lOT and F = 2 x 10~^. Note that the azimuthal angle (j) iii 
Fig. 6 is shifted by 45° compared to Fig.l because of the difference of the unit cells. (In our first principles 
calculations, the unit cell contains two Fe atoms and two As atoms. Meanwhile in the five band model of 
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ref.23, the unit cell contains one Fe atom. The Fe-Fe bond direction differs by 45° between the two unit 
cells.) If we compare Fig. 6(b) with Fig. 1(b), they are qualitatively the same. The peaks at = 45° and 
(j) = 135° in Fig.l(b) are observed at = 0°, 90°, and 180° in Fig.6(b). The peak at (/> = 90° in Fig.l(b) 
is observed at = 45° and = 135° in Fig. 6(b). However, the peak widths are different. The discrepancy 
is mainly associated with the Fermi surface shape difference between FS4 and FS5 in Fig. 1(a) and FS4 
and FS5 in Fig. 5(b). In addition there arc non-negligible discrepancies in the Fermi velocities between the 
five band model and the first principles calculation. Those discrepancies lead to the differences of each 
band's weight in the interlayer conductivity. The Fermi velocities are the derivatives of the dispersion 
energies with respect to kx or ky. However, we may not expect that the five band model should reproduce 
derivatives of the energy dispersions of the first principles calculation. Although there is some quantitative 
differences, the five band model reproduces the first principles result qualitatively. 

In hole doped case, the result of the interlayer conductivity is similar to that of n = 6: For 10% hole 
doping case (n = 5.9) the Fermi surface Fig. 5(a) is similar to n = 6 case (Fig. 5(b)) and so their interlayer 
conductivities are similar (Figs. 6(a) and 6(b)). Upon further hole doping, the electron Fermi surfaces FS4 
and FS5 shrink. But the shrink of those Fermi surfaces does not lead to qualitative difference of the 
interlayer conductivity. 

In the electron doped side, for n = 6.2, the hole Fermi surface around (vr, vr) (FS3) disappears (Fig. 5(c)). 
The effect of this disappearance of the Fermi surface FS3 is clearly seen in the interlayer conductivity 
Fig. 6(c): The peak at <^ = 45° and (j) = 135° in Figs. 6(a) and 6(b), which is associated with the flat regions 
of the Fermi surface FS3, vanishes in Fig. 6(c). For n = 6.4 the two hole Fermi surfaces around (0,0) (FSl 
and FS2) disappear, and the peak at (/> = 90° in the interlayer conductivity becomes sharp (Fig. 6(d)) 
because of large flat regions of FS4 and FS5 (Fig.5(d)). 

5. Summary and Discussion 

In this paper, we have proposed a interlayer conductivity formula under an in-plane magnetic field 
to study the Fermi surface topology of layered multiband systems. Applications of the formula to one 
of the FeAs-based superconductors, LaFeAsOi-jjF-c, have been discussed. Although the total interlayer 
conductivity is given by the sum of all band contributions, it is possible to extract information about 
the Fermi surface topology from the peak positions in the interlayer conductivity as a function of the 
azimuthal angle of the magnetic field. Flat regions of the Fermi surface give rise to peaks in the interlayer 
conductivity, while nearly circular Fermi surfaces do not show oscillating patterns. We have discussed the 
difference between the paramagnetic state and the magnetic state of the parent compound LaFeAsO. As 
for the presence of the magnetic state, the oscillation period is 180° for the magnetic state but 90° for 
paramagnetic states. To study the doping effect, we have taken the five band model of ref.23 and discussed 
the Fermi surface topology change upon doping and how such changes are reflected to the interlayer 
conductivity. 

In the analysis above, the interlayer coupling has been taken into account through the interlayer current. 
The energy dispersion along the z axis is not included. If the electron band dispersion along the z axis is 
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Fig. 5. (Color online) Electron density dependence of the Fermi surface calculated from the five-band model by 
Kuroki et al?"^ 

simply represented by an additional term of e^*^^ = —2tcCOskz, then the Lorentz function factor in eq.(lO) 
has dominant contributions at = 0, vr. However, if the Fermi energy is much larger than t^., the Fermi 
surface topology change is negligible. For such a simple energy dispersion form in the /c^-direction, one may 
apply the Shockley's tube- integral formula^^ to analyze the angle-dependent magnetoresistance oscillations. 
The angle-dependent magnetoresistance peaks are connected with the vanishing of the electronic group 
velocity perpendicular to the layers.^'' Since the vanishing condition depends on the Fermi wave-vector, 
which is azimuthal angle dependent, it is possible to extract further information about the Fermi surface 
topology. ^'^ 

As for the scattering, we have assumed the same scattering time for all bands for simplicity. However, 
the scattering time could be different for each band. In this sense, the measurement of the interlayer 
magnetoresistance is useful for the determination of the scattering time, because oscillating patterns are 
discernible from a certain magnetic field. Although the onset is not directly related to F, one can extract 
r by comparing first principles calculation results with the experiments. 

The results above have been obtained at zero temperature. In order to evaluate the finite temperature 
effect precisely, we need to sum over contributions from energy contours with including the first derivative 



8/12 



J. Phys. Soc. Jpn. 



Full Paper 




20 40 60 80 100 120 140 160 180 
iMdefl) 



(a) n = 5.9 




0. 

0. 
0. 
0. 
0. 
0. 
0. 
0. 

0.1 





20 40 60 SO 100 120 140 160 ISO 



1.1 
1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 








Total 




FS1 




FS2 




FS3 




FS4 ; 




FS6 





20 40 60 SO 100 120 140 160 18 
iKdeg) 



(b) n = 6 



— 0.6 

"a 




20 40 60 80 100 120 140 160 18 
(Kdeg) 



(c) n = 6.2 



(d) n: 



Fig. 6. (Color online) Electron density dependence of the Fermi surface calculated from the five-band model by 
Kuroki et al?"^ 



of the Fermi-Dirac distribution function with respect to the energy. The finite temperature effect can be 
estimated for a simple energy dispersion case. If the energy dispersion is well approximated by = h'^{k'^ + 
Xky)/{2m), then the amplitude of the conductivity is somewhat suppressed for finite temperatures but the 
characteristic temperature for the suppression is scaled by the Fermi energy. Furthermore, the temperature 
dependent factor is factorized out from the magnetic field dependent part. For energy dispersions with 
complicated kx and ky dependences, this simple argument does not apply. However, we expect that the 
finite temperature effect on P is more significant than that on the other contributions in the interlayer 
conductivity. 

The interlayer magnetoresistance measurement can be used to find nodes in the superconducting state 
if there are nodes in the gap as discussed by Bulaevskii et al?^ If the z axis current exceeds the Josephson 
critical current, then one can measure the interlayer magnetoresistance. In the presence of the gap nodes, 
the interlayer magnetoresitance detects the nodes if the applied magnetic field is parallel to the normal 
state Fermi velocity at the node. However, the formula given in ref.28 is not directly applicable to the 
sib-wave gap proposed for FeAs-based superconductivity.^^' 

Finally we comment on advantages of using the interlayer magnetoresistance for investigating Fermi 
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surface topology. Although ARPES is a powerful measurement to extract information about Fermi surface 
topology, the interlayer magnetoresistance experiments are applicable to more compounds than ARPES, 
since it is not necessary to cleave a sample to obtain a clean surface. In addition, the interlayer magnetore- 
sistance is associated with the bulk properties. Therefore, the measurement is not sensitive to the surface 
condition. There are now many FeAs-based superconductors. In order to make clear common underlying 
electronic properties for the mechanism of superconductivity, we need to compare two different systems. 
In this sense, the interlayer magnetoresistance is useful because of its wide applicability. 

We would like to thank T. Shibauchi for drawing our attention to ref.28. This work was supported 
by the Grant-in- Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and 
Technology (MEXT) of Japan, the Global COE Program "The Next Generation of Physics, Spun from 
Universality and Emergence," and Yukawa International Program for Quark-Hadron Sciences at YITP. 
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